Nonuniversal routes to universality: Critical phenomena in colloidal dispersions 
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We investigate critical phenomena in colloids by means of the renormalization-group based hier- 
archical reference theory of fluids (HRT). We focus on three experimentally relevant model systems, 
namely the Asakura-Oosawa model of a colloidal dispersion under the influence of polymer-induced 
attractive depletion forces; fluids with competing short-range attractive and longer-range repulsive 
interactions; solutions of star-polymers whose pair potential presents both an attractive well and an 
ultrasoft repulsion at shorter distance. Our results show that the ability to tune the effective inter- 
actions between colloidal particles allows one to generate a variety of crossovers to the asymptotic 
critical behavior, which are not observed in atomic fluids. 
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The description of colloidal dispersions in terms of 
effective interactions acting between the macroparti- 
cles [H has opened the realm of complex fluids to mod- 
ern Liquid-State Physics, and is the main reason for the 
rekindled interest in this discipline. On the one hand, 
using colloidal particles as macroatoms allows visual in- 
spection of the fluid, sohd and glassy states of matter 
which occur also in atomic substances 0] . On the other 
hand, the diversity of effective interactions, as well as 
the possibility of tuning them by changing the parame- 
ters which characterize the dispersion, leads to a phase 
behavior which is much richer than that encountered in 
atomic fluids 0. This ability of engineering the phase 
diagram Q is fascinating both in its own sake, and for 
its potential in technology. Compared to the study of 
their overall phase behavior, critical phenomena in col- 
loids have received less attention. A likely reason for this 
situation is that, in contrast with the menagerie of pos- 
sible phase equilibria, the most relevant feature of criti- 
cal phenomena is their universality, as expressed by the 
critical indexes being the same for all the systems be- 
longing to the same universality class. For fluid-fluid 
(FF) transitions, this is expected to be the Ising one for 
both simple and colloidal fluids such as protein solutions 
and colloid-polymer mixtures, so it might be argued that 
there is little to learn in going from simple fluids to col- 
loids. However, while universality comes along asymp- 
totically close to the critical point, the approach to the 
asymptotic regime is by no means universal: how the 
asymptotic critical regime is reached, and how far from 
the critical point critical fluctuations are expected to be- 
come important, are relevant questions whose answers 
can be quite different from those found in atomic fluids, 
since the specific features of colloid-colloid effective po- 
tentials may induce peculiar crossovers, and affect the 
size of the critical region. 

In this Letter, we show how the onset of critical fluc- 



tuations is affected by the specific form of the interac- 
tion by focusing on three model systems, all of which 
have received much attention because of their phase be- 
havior and structural properties: the Asakura-Oosawa 
(AO) pair potential for hard-core colloidal particles in the 
presence of non-adsorbing polymer in solution; a hard- 
core two- Yukawa (HCTY) fluid with competing short- 
range attractive and longer-range repulsive interactions; 
a model of a star-polymer (SP) solution, where the ul- 
trasoft repulsion between the stars is followed by an at- 
tractive well. We have investigated their critical behavior 
by the hierarchical reference theory (HRT) In HRT, 
the attractive part of the interaction in Fourier space is 
gradually switched on by introducing an evolving infrared 
cut-off Q such that at any intermediate stage of the pro- 
cess, fluctuations with characteristic length L > Q are 
strongly suppressed. The fully interacting system and 
the onset of long-range correlations are then recovered 
in the limit Q ^ 0. HRT is optimally suited to the 
present study because it implements the renormalization 
group (RG) method, and therefore yields a realistic, non 
mean-field-like treatment of the long- wavelength fluctua- 
tions which drive criticality and phase separation. On the 
other hand, this approach keeps all the information about 
the non-universal properties which are expected to de- 
pend on the speciflc features of the interaction, and would 
be lost in a coarse-graining approach. The crossover to 
the asymptotic power-law behavior can be described by 
introducing an "effective" critical exponent, deflned as 
the local slope of the quantity in hand as a function of 
the reduced temperature i in a log- log plot. In the follow- 
ing, we will adopt the standard dcflnition t = {T—Tc)/Tc, 
Tc being the critical temperature. 

We first consider the AO model. The AO effective pair 
potential of polymer-induced depletion forces Q consists 
of an attractive tail which vanishes identically beyond the 
polymer diameter, measured from the hard-core surface. 



2 



1 

0.8 

0.7 
1 

> 
> 

0.8 



• AO, q=0.4 
aAO, q=0.56 
□ AO, q=0.8 
X Yukawa, z=1.8 



□ A 
□A 



X X ° O^n'^ ^ AA A A 



-6 -5 -4 -3 -2 -1 
log,„(t) 

FIG. 1: AO effective exponents 7efi (upper panel) and v^s 
(lower panel) for Xred and 5, normalized by their asymptotic 
values, compared with the effective exponents of an attractive 
HCY potential with inverse-range parameter z = 1.8. For the 
AO potential, the reduced effective temperature t is {rfp^^^^t ~ 
rfp)/rfp, where rfp is the reservoir packing fraction. 



The strength of the interaction is controlled by the fugac- 
ity Zp of the polymer or, equivalently, by the packing frac- 
tion of the pure polymer solution in osmotic equilib- 
rium with the colloid-polymer mixture. This "reservoir" 
packing fraction plays the role of the inverse temperature 
in a thermal fluid. Consistently with the above definition 
oft, we have set t = {Vp,crit -Vp)/'np, ^p.crit being the crit- 
ical value of 77p. A HRT study of the overall phase dia- 
gram and asymptotic critical behavior of several thermo- 
dynamic quantities has been presented elsewhere 0, [1] ■ 
However, the approach to criticality has not been studied 
before. In Fig. [1] we have plotted, for several polymer- 
to-colloid size ratios q, the effective exponents ^eS, i^eff 
normalized by their asymptotic values 7, z^, which de- 
scribe the divergence of the (reduced) compressibility 
Xrod and correlation length ^ on the critical isochore as Tc 
is reached. In HRT one has 7 = 1.378, ly = -//2 = 0.689. 
Figure [1] also shows the effective exponents of a hard-core 
Yukawa (HCY) potential with attractive tail and inverse- 
range parameter z = 1.8, which is appropriate to model 
the pair potential of simple atomic fluids (inverse ranges 
are in units of cr^^, a being the hard-sphere diameter). 
As expected on the basis of the RG predictions [§], as 
q is decreased, the asymptotic critical regime is reached 
at larger reduced temperatures, and for q = 0.4, 7cff is 
about 95% of its asymptotic value already at reduced 
effective temperatures t ~ 0.1. This is different from 
what is observed in the HCY fluid, whose 7off at these 
reduced temperatures is still mean-field-like. We did not 
consider size ratios q < 0.4, since our implementation of 
HRT is not suited to very narrow potentials. Moreover, 
in the AO potential, FF phase separation is preempted 
by freezing for q < 0.4 0, [i31j even though at smaller q it 



can still be observed as a metastable transition. For the 
values of q considered here, the AO depletion interaction 
has many-body contributions, which arc not accounted 
for by the effective pair potential which we have used. 
However, the critical behavior is captured satisfactorily 
already at the pair potential level [Sj]. On the basis of 
our results, we expect AO colloid-polymer mixtures to 
exhibit nearly asymptotic Ising scaling at reduced effec- 
tive temperatures t ~ 0.1. We also note that there is no 
exponent renormalization due to the binary nature of the 
system, provided the reservoir packing fraction is used as 
the relevant field driving criticality. 

In a recent paper [lH , the critical behavior of the corre- 
lation length of colloid-polymer mixtures has been inves- 
tigated experimentally, and the critical exponent v was 
found to be close to the mean-field value ly = 0.5. On the 
other hand, Monte Carlo studies supplemented by finite- 
size scaling unambiguously proved that the AO model be- 
longs to the Ising universality class [1, [l3] ■ This rises 
an interesting question about the possible presence of an 
unforeseen long-range interaction in the mixtures studied 
in Ref. [ll|. 

Besides the tunable short-range attractions of which 
the AO potential represents the best-known example, flu- 
ids where the short-range attraction is followed by a re- 
pulsion which takes over at long distance have been re- 
ceiving increasing attention. It has long been acknowl- 
edged [l3| that the competition between attraction and 
repulsion leads, for strong enough repulsion, to the dis- 
appearance of bulk FF phase separation. As the tem- 
perature is lowered, this is replaced by the occurrence 
of nonhomogeneous phases that may have different mor- 
phologies [l5|, e.g. cluster- or stripe-shaped. 

Instead of considering the nonhomogeneous phases, 
here we focus on a regime of relatively weak repulsion, 
such that the competition is not so strong as to cause 
the disappearance of FF phase separation, and investi- 
gate the effect on the critical behavior. Our model in- 
teraction is the HCTY potential with inverse-range pa- 
rameters zi = 1, Z2 = 0.5 for the attractive and repulsive 
contribution respectively [l^l ■ The relative amplitudes A 
of the repulsion is close to the value above which the FF 
transition disappears, and microphase formation is ex- 
pected. This regime is characterized by an extremely 
flat FF coexistence curve [l^. In Fig. [21 Xrcd and ^ 
have been plotted together with the effective critical ex- 
ponents. The most striking effect is the steep increase of 
Xred and ^ at reduced temperatures t ~ 0.1, with a corre- 
sponding sharp crossover of the effective exponents. The 
enhancement of large density fluctuations can be traced 
back to the tendency towards cluster formation triggered 
by the competition. A sign of this "incipient cluster- 
ing" is that, unlike in Lennard- Jones (LJ)-like fluids with 
purely attractive tail potentials, one flnds damped oscil- 
latory decay for the correlations even for states very close 
to critical 17[. On further approaching the critical point. 
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FIG. 2: Reduced compressibility (upper panel) and corre- 
lation length (lower panel) in units of the particle diame- 
ter of a HCTY fluid with competing short-range attraction 
and longer-range repulsion. The insets show 7eff, VeS com- 
pared with the effective exponents of a HCY potential {ucS 
for A = 0.096 has been omitted for clarity). See text for the 
meaning of the lengths A and ^. In the oscillatory regime, 
only the real part of A and ^ has been shown. Lines are a 
guide for the eye. 



the decay becomes monotonic, but takes a double expo- 
nential form described by two characteristic lengths A, 
fi l3j A similar behavior has in fact been found in 
many frustrated systems [l9j . In our treatment, A and fi 



are obtained from expanding the reciprocal of the struc- 
ture factor S{k) at small wave vector k up to fc^. Asymp- 
totically close to the critical temperature, A reduces to 
the usual Ornstein-Zernike (OZ) correlation length, while 
fi saturates at a finite value. The change from oscilla- 
tory to monotonic behavior is marked by the bifurcation 
in Fig. [21 where A, fi turn from complex conjugate to 
real. This is an example of the well-known Kirkwood 
crossover [2^. We observe that j/cff shows both a peak 
at a temperature slightly higher than that of the bifur- 
cation, and a divergence at the low-temperature side of 
the bifurcation, where A takes a vertical slope. The lat- 
ter feature is due to the coalescence of A and /i at the 
Kirkwood line. 

It is tempting to intepret the above scenario as 
a cluster-driven criticality, whereby particles are first 
strongly correlated within regions of size ~ /i, which in 
turn become correlated over length scales of size ~ A 
as the critical point is reached. In the asymptotic criti- 
cal regime, Xred diverges as A^ like in the OZ picture, but 
with an amplitude ~ /i^, much bigger than that of LJ-like 
fluids, see Fig.O The possibility of clusters behaving like 
super-atoms has already been pointed out in the context 
of the glass transition [2l|. In our case, however, clusters 



log,o(t) 

FIG. 3: Same as Fig. [l]for three SP solutions. Lines are a 
guide for the eye. 



denote regions of enhanced particle-particle correlations, 
rather than stable particle aggregates. 

The anomalous crossover of the HCTY model close 
to the stability limit of the liquid-vapor transition has 
not been observed experimentally yet. However, it seems 
likely that such a regime can be experimentally reached. 
As an example, let us consider the system studied by 
Cardinaux et al. [2^ as a model of lysozyme solutions. 
They investigated cluster formation for volume fractions 
Tj of lysozyme in the range 0.085 < 77 < 0.201 in a so- 
lution of sodium hydroxyde at concentration c = 8mM. 
An estimate of the electrolyte concentration below which 
the liquid-liquid transition is expected to disappear can 
be obtained within the random phase approximation by 
locating the concentration at which the convexity of the 
Fourier transform of the potential at zero wave vector k 
turns from positive to negative. According to this crite- 
rion, c = BmM is indeed well inside the cluster region, 
while at ij ~ 0.201, the threshold concentration is at 
c ~ 56mM. 

Star-polymer solutions differ from the two systems con- 
sidered above inasmuch as the interaction which accounts 
for the excluded volume effect is not the hard-sphere po- 
tential, but an ultrasoft repulsion of entropic origin with 
a logarithmic core and a repulsive Yukawa tail [23| . The 
size of the star a is of the order of the radius of gyra- 
tion. It has been shown that, when an additional attrac- 
tive interaction between the stars is present, the solution 
may display two FF transitions, each ending at a crit- 
ical point [2J|. Here we focus on the behavior in the 
neighborhood of the lower-density critical point, which 
we expect to be more easily accessible by experiments. 
We considered three values of the arm number or "func- 
tionality" /, namely / = 12, / = 24, / = 32. The at- 
tractive well was modeled by the same Fermi function 
used in Ref. [3] . Figure [3] shows the effective exponents 
7eff, t'eff- The main qualitative difference with respect to 
the simple-fluid case represented by the HCY fluid with 
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z = 1.8 consists in the non-monotonic behavior of the 
effective exponents. As one moves away from the crit- 
ical temperature, v^,^ shows a minimum in the interval 
10~^ < t < 10~^, and increases sharply as t is further 
increased. This is due to the effect of the soft repulsion 
which, at the small densities considered here, is more sub- 
stantial than in a hard-core system of comparable size. 
Moreover, both 7eff and show another minimum at 
t ~ 10^^ and bend towards their asymptotic values as t is 
further reduced. The accuracy limit entailed by our nu- 
merical calculation prevented us from reaching reduced 
temperatures smaller than t ^ 10^®, so that the asymp- 
totic exponents could not be observed. Notwithstanding 
this, the quantitative deviations from simple-fluid behav- 
ior for t < 10~^ arc rather small, and the effective ex- 
ponents remain well apart from their mean-field values 
even in the region of the minimum, at least for the val- 
ues of / investigated here. This behavior is different from 
what is found in the linear-polymer solutions considered 
in Ref. [2^. In that case, the correlation length ^ away 
from the critical point is considerably smaller than the 
polymer size a. As the critical point is approached and 
^ becomes comparable to a, the solution shows a rather 
sharp crossover from mean-field to Ising critical behav- 
ior. Here, on the other hand, the size of the star a is the 
natural lengthscale of both the repulsive and attractive 
interaction, so that ^ is of the order a even away from the 
critical point. This requires that the quality of the sol- 
vent should be unchanged at the transition, so that the 
inter-star attraction is not overwhelmed by the attrac- 
tion between monomers within the same star. Experi- 
mentally this condition can be reproduced by considering 
effective depletion attractions between stars in multicom- 
poncnt star-chain [26| or star-star mixtures pTj . These 
systems have been looked upon for a big component in 
solution with higher / than those we investigated here. 
Nevertheless, from the trends evidenced and the compar- 
ison between the effective interaction shapes, we expect 
the appearence of a stable low-density FF transition for 
highly or moderately asymmetric mixtures with small / 
(/ ^ 50) for the big component. 

In summary, we have presented an investigation of the 
critical behavior of three model fluids, whose pair inter- 
actions arc relevant for a wide class of complex liquids, 
although they certainly do not aim at covering the whole 
subject of criticality in colloidal systems. In the AO 
model, nearly asymptotic Ising-likc behavior is observed 
in a wide interval of reduced effective temperatures for 
the polymer-colloid size ratios at which the effective AO 
pair interaction is a reliable representation of the binary 
colloid-polymer mixture. In fluids with competing inter- 
actions, the enhancement of density fluctuations due to 
competition causes strong crossovers with anomalously 
high effective critical exponents j^s, t'off and critical am- 
plitudes, as well as the emergence of two characteristic 
lengths. In star-polymer solutions with star-star attrac- 



tion, we find a non-monotonic behavior of the effective 
exponents 7eff, VeS close to the critical temperature. We 
expect that these predictions can be tested by experi- 
ments in colloid-polymer and polymer-polymer mixtures 
and protein solutions. 
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